Effects of Rising Angle on Upstream Blades and Intermediate Turbine Duct 
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With the improvement of requirement, design and manufacture technology, aero-engines for the future are cha- 
racterized by further reduction in fuel consumption, cost, but increment in propulsion efficiency, which leads to 
ultra-high bypass ratio. The intermediate turbine duct (ITD), which connects the high pressure turbine (HPT) with 
the low pressure turbine (LPT), has a critical impact on the overall performances of such future engines. There- 
fore, it becomes more and more urgent to master the design technique of aggressive, even super-aggressive ITDs. 
Over the last years, a lot of research works about the flow mechanism in the diffuser ducts were carried out. Many 
achievements were reported, but further investigation should be performed. With the aid of numerical method, 
this paper focuses on the change of performance and flow field of ITD, as well as nearby turbines, brought by 
rising angle (RA). Eight ITDs with the same area ratio and length, but different RAs ranges from 8 degrees to 45 
degrees, are compared. 

According to the investigation, flow field, especially outlet Ma of swirl blade is influenced by RA under po- 
tential effect, which is advisable for designers to modify HPT rotor blades after changing ITD. In addition to that, 
low velocity area moves towards upstream until the first bend as RA increases, while pressure loss distribution at 
S2 stream surface shows that hub boundary layer is more sensitive to RA, and casing layer keeps almost constant. 


On the other hand, the overall total pressure loss could keep nearly equivalent among different RA cases, which 


implies the importance of optimization. 
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Introduction 


Efficiency and SFC(Specific Fuel Consumption) are 
two important indicators of engine. In order to raise the 
performance, and to decrease the cost, higher efficiency 
and lower SFC are goals which have been pursued unre- 
mittingly. To fulfill such requirement, larger bypass ratio 
engine received more attention for its advantage and po- 
tential. However, when bypass ratio increased and hence 
the power needed by the fan remains, it is logical that the 
radius of LPT blades is tend to be larger. So the radius 
difference between HPT and LPT increased on certain 


extent. However, such radius difference might introduce 
some worry from the designer’s point of view. ITD, so, 
not only plays the role of connecting the HPT and the 
LPT, but also with the function of diffusing air to de- 
crease Ma with lower pressure loss, is believed to be in- 
vestigated essentially in more detail. 

The performance of ITD can be influenced by many 
factors, for example, area ratio(AR), non-dimensional 
length(ratio of ITD length to inlet height), and so on. 
Among them, upstream flow condition is an apparent 
factor. Dominy et al.!"! investigated the influence of swirl 
upon the flow through a diffusing duct. On the same ex- 


Nomenclature 

AR Area Ratio Abbreviation 

HPT High Pressure Turbine Cps Static Pressure Coefficient 
ITD Intermediate Turbine Duct Cpt Total Pressure Coefficient 
LPT Low Pressure Turbine Pso Inlet Static Pressure 


Ma Mach Number 
RA Rising Angle 

SP Static Pressure 
TP Total Pressure 


perimental facility, Dominy et al.”! investigated the in- 


fluence of the upstream blade wakes on the performance 
of the duct. To the experience of the author, their work is 
the first experiment dealing with the pseudo-steady in- 
fluence of two-dimensional wakes extending from hub to 
the casing on ITDs. More recently, Hu et al.”! and Zhang 
et al.'*! researched the influence of casing swirl and hub 
swirl distributions separately on the flow field within 
ITDs with experimental and numerical methods syste- 
matically. Miller et al.'°! published their work of re- 
searching the migration and dissipation of flow pheno- 
mena in a swan-necked diffuser downstream of a tran- 
sonic high-pressure turbine stage. And in the next year, 
Miller et al.!°! extended the knowledge of this area, that 
they investigated the interaction mechanisms between a 
high-pressure rotor and a downstream vane located in a 
diffuser. The influence of upstream HPT rotor gap size 
variation on the ITD’s flow field is discussed by Marn et 
al.!”!°l Göttlich et al.!, and Sanz et al.”!. In addition to 
that, Axelsson et al.!!"! revealed the complex flow struc- 
ture development within ITD. 

Gottlich et al.!'7! and Pullan et al. "have studied the 
flow field within ITD. Johansson et al.''*! finished a work 
that focused on a joint experimental and numerical inves- 
tigation of an intermediate duct configuration. Marn et 
al.l"! compared the flow through two different s-shaped 
intermediate turbine duct configurations with different 
diffusion ratio. Hu et al.''®! optimized a baseline ITD 
geometry by the software of Numeca to get a second ITD 
with 25% larger mean RA and a third ITD with 20% 
larger area ratio. 

Gottlich!'” summarized the previous work about ITD. 
The author reviewed the flow evolution through interme- 
diate turbine diffusers and discussed the influence of dif- 
ferent effects in a systematic way. This was a meaningful 
work for all the following investigations to have a good 
start and get a correct direction when they try to develop 
ITD. In addition to that, other useful investigations!'*!?*! 
about the ITD were finished by researchers. 

Upstream flow condition is one important and appar- 
ent factor for ITD performance, which could also be rea- 
lized by means of experimental facility and receives 
much study as above introductions. Another one, which 
is vital as well and attracts much attention but less inves- 
tigation to be done, is the duct geometry itself. 

Norris and Dominy”” used both experimental and nu- 
merical results to show the differences between two S- 
shaped ducts with the same area ratio and radial offset 


Ps Local Static Pressure 
Pto Inlet Total Pressure 
Pt Local Total Pressure 


but with different axial length. 

Couey et al.”®! developed the original annular-dif- 
fuser performance chart for ITDs This original chart has 
been widely used for diffuse-type duct design after So- 
vran and Klomp'”*! reported it. With numerical method, 
Couey et al. assessed the interaction of duct slope and 
pitch-wise turning with area ratio and length using ap- 
proximately sixty different duct configurations. After 
investigation, they got the conclusion that the influence 
of slope and exit pitch angle can be as significant as the 
area ratio and length parameters traditionally used to 
correlate duct performance. And a combination of new 
area ratio-length and slope-length parameters was found 
to segregate ducts between separated and non-separated 
cases. At last, they modified the Sovran and Klomp 
Chart. 

A more detailed investigation of the effects of mean 
RA and area ratio on the aerodynamics of ITDs was done 
by Zhang!*”! et al. in the large scale, low-speed, annular 
wind-tunnel facility at the National Research Council of 
Canada. After study, they found that the duct mean rise 
angle determined the severity of adverse pressure gra- 
dient in the casing’s first bend whereas the duct area ratio 
mainly governed the second bend’s static pressure rise. 

In above contents, Norris and Dominy did the work to 
research the flow field of two ducts with the same area 
ratio and radial offset, but these ducts had different axial 
lengths. A more detailed work was done by Couey et al.. 
Almost sixty models were simulated to find the influence 
of duct slope parameter and pitch-wise turning parameter 
on ITD, they focused on the relationship between overall 
performance and these key design parameters of ITD, but 
without further investigation of the flow field detail cor- 
respondingly. In addition to that, as author pointed out, 
without upstream blades, the real coupled influence was 
not captured. Because upstream wake is the main source 
of pressure loss, that should be paid more attention to if 
one wants to evaluate the ITD performance more accu- 
rately. Another question confused the writer is that 
whether the throat area of LPT vane keeps unchanged, 
which should be a more important parameter for vanes. 

What’s more, almost all the attentions were focused 
on the performance of ITD, while how the nearby tur- 
bines, especially upstream blades, were influenced re- 
ceived much less focus. 

In this paper, ITD model is designed with upstream 
swirl vanes (which can supply wake and inlet swirl an- 
gle), and downstream LP nozzle vanes (which can supply 


real downstream flow field with effect of potential and 
blockage), as Figure 1 shown. Eight configurations are 
chosen with the aid of numerical method. All of these 
model configurations have the same AR and axial length, 
but different radial offsets. Performance comparison and 
flow field analysis for these models will be presented in 
the current paper with the purpose of identifying the in- 
fluence of the radial offset of ITD on turbine perfor- 
mance and flow field. It is hopeful that the result would 
give some sensible thoughts when more aggressive ITD 
design is ready to be used in the future advance engine. 


Duct Inlet 


Fig. 1 Defination of ITD model 
Numerical Model and Boundary Conditions 


The original conservative ITD prototype is a well- 
designed duct which is positioned between HPT and LPT 
with RA of 8 degrees from an aero-engine. With the pur- 
pose of further improving power output of LPT, it is pre- 
pared to study the ultimate potential of such ITD in this 
paper. 

In order to investigate the influence of RA on ITD 
flow field separately, apart from the RA, other parame- 
ters are kept unchanged in these models. Inlet conditions 
and swirl blades are kept the same, as well as the axial 
outlet position (inlet of LPT nozzle). The parameters of 
ITD and LP nozzle are listed in Table 1. The numerical 
models with RA of 8(A8), 15(A15), 20(A20), 25(A25), 
30(A30), 35(A35), 40(A40) and 45(A45) degrees are 
investigated in parallel. However, the inner and outer 
walls of ITD’s outlet section are modified to make sure 
that the AR of ITD is equivalent. Meanwhile, the hub and 
casing curvature of LPT nozzle are adapted to the mod- 
ified duct. As Figure 2 shown, A8, A20, A30 and A40 are 


presented. 

Table 1 Parameters of ITD 
AR L/H Inlet Swirl Angle(°) Throat Area of LPT(m’) 
1.57 3.38 20 0.0198 


Fig.2 The sketch of different RA ITDs 


What should be paid more attention is that the LPT 
nozzle blade numbers of all the models are kept the same 
(44 for swirl blades and LPT nozzle). As what has been 
mentioned above, in the current paper, authors would like 
to distinguish the RA’s influence on the ITD. To fulfill 
such goal, it is believed by the authors that the same 
throat area of the LPT nozzles must be realized to guar- 
antee the equivalent mass flow rates. Usually there are 
two ways to realize such requirement. One way is in- 
creasing the blade number of nozzle with all the blades 
profile kept unchanged; another way is altering the blade 
profile on certain extent, on one section or multi sections, 
while the blade number keeps unchanged. Although it is 
easier to finish the model modification with the first me- 
thod, the unexpected relative larger weight of LPT nozzle 
accompanied with the increase of profile loss may elimi- 
nate our maneuver for further improving the turbine per- 
formance. So the second strategy is adopted. The throat 
area remains constant of 0.0198(m”) in all the studied 
cases. And this is achieved by altering the nozzle’s outlet 
mental angle, stagger angle and ratio of outlet passage 
width to the thickness of blade. Unguided turning angle 
is amended to smooth the blade surface. The boundary 
conditions for all models are the same, with a constant 
in-mass-flow of 8.52 kg/s at inlet, and back pressure of 
240 kPa. So the Reynolds number is kept 310° to avoid 
its influence on boundary layers. 

Commercial CFD software of Numeca is used. The 
topology of building model mesh is O4H, and matching 
periodicity technology is used to get an accurate flow 
field in blade passages. Figure 3 is the numerical model 
and structured grids of A40 as an example created by 
IGG. The mesh at blade root and leading edge is magni- 
fied locally to show the detail mesh. The orthogonality of 
mesh is good enough to reflect the flow field. And Y+ 
near wall is set below 1 to fulfill the requirement of the 
Shear Stress Transport(SST) turbulence model. With the 
same count of blades, full non matching frozen rotor be- 
tween swirl vanes and LPT nozzle is chosen to provide 
steady upstream wake for ITD. 

The quantity of model mesh is an important factor for 
numerical results. So a test of mesh independence to si- 
mulation result is necessary. Figure 4 is the relationship 
between total pressure ratio and mesh quantity of A8. 
This regularity applies to other models as well. From this 
chart, the mesh count will not influence the pressure ratio 
when it is more than 3 million, which means that the flow 
field is independent of mesh count. So each model will 
be simulated with more than 3 million meshes and their 
computation results are ready to be compared in the cur- 
rent paper. 
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Fig.3 Numerical Model and Mesh of A40 
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Fig. 4 Test of Mesh Independence 
Results and Discussion 


All the models are simulated while RA is changed. In 
the following contents, the numerical results will be ex- 
hibited systematically. According to the simulation result, 
RA is an important factor which can change the perfor- 
mance of ITD and nearby turbine blades. 

After simulation, separation occurs in A45, and could 


not be diminished by modifying LPT nozzle blade profile. 


So it is advised in this paper that RA should be set small- 
er than 45 to prevent separation. A45 will not be dis- 
cussed to avoid severe numerical error. 

Figure 5(a) and Figure 5(b) are yaw angle and Ma at 
ITD’s inlet section, Plane Cl. 
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ITD Inlet Angle 


Passage Height 
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(a) ITD Inlet Yaw Angle 
ITD Inlet Ma 


(b) ITD Inlet Ma 


Fig. 5 


From Figure 5(a) and Figure 5(b), compared with Ma, 
yaw angle is less influenced accompanied with the varia- 
tion of RA, which means that angle is more stable under 
the influence of potential field. Except near hub region, 
the initial A8 case has a larger yaw angle in the main 
flow region than other cases. This can be deduced by the 
fact that the relatively smaller concave curvature along 
inner wall of A8 brings such phenomenon. In A8, more 
specifically, the tangential velocity is smaller near hub 
region, which reduces the yaw angle slightly. As the same 
reason, yaw angle above 50% passage height of A40 
model is the smallest, especially at 95% passage height. 

At the ITD inlet section it can be evidently found that 
larger RA accompanied with smaller Ma. While A40 case 
has evident Ma unevenness, its counterpart A8 case 
owned more even Ma model. So it is important to focus 
on this characteristic to get a more even flow field. 

The detailed flow field in S3 stream surface is useful 
to observe the local flow characteristic in totality. Figure 
6 presents the ITD inlet (Plane Cl) Ma. Within the 
smaller RA case, high Ma region almost occupies all the 
passage height. However, the range of low Ma increases 
gradually when RA increased steadily. As RA reaches its 
peak in case A40, Ma, higher than 0.7, can be only seen 
above about 85% passage height. Stemming from swirl 
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blade’s pressure side, low velocity region developed and 
the core of low energy fluid remained closed to its pres- 
sure side. 

According to Ma and yaw angle variations, ITD RA 
affects HPT rotor deviation angle. As the result, HPT 
performance is influenced, which should be forecasted 
during turbine design process. 

The SP distributions on swirl blade surfaces are com- 
pared in Figure 7, and black dotted lines represent SP 
isolines. According to Figure 7(a), the blade surface SP 
increases when RA is larger. High pressure(larger than 
292kPa) regions extend from area near casing to entire 
blade surface. So the fluid decelerates near blade in larg- 
er RA model. Low pressure corner(less than 254kPa) 
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Fig. 6 ITD Inlet Ma (From Downstream) 
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Fig. 7 SP Distribution on Swirl Blade Surface 
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regions locate in blade trailing edge near root in A8. 
However, it moves near casing in A30 and A40. As the 
swirl blade passage is horizontal, and blade profile was 
designed without twisting, the pressure distribution is 
influenced by stream-wise adverse pressure gradient. 
Low pressure corner area will change blade deviation 
angle. The feature of suction surface pressure is various. 
Low pressure corner region is not observed instead of 
high pressure corner region developing from area near 
casing of A8 to area near root of A40. Low pressure re- 
gions narrow down from entire blade surface of A8 to 
trailing edge area of A40 (the blank area at trailing edge). 

In subsonic turbine stage, on one hand, more attention 
should be received for downstream flow field; On the 
other hand, different RA can cause change to upstream 
flow field simultaneously. With larger ITD outlet radius 
of investigated models, the inlet flow for LPT nozzle 
differs. The yaw angle and velocity distribution would 
present various features to affect downstream turbine 
performance. 

More specifically, flow field in S2 surface of duct will 
be discussed. If separation occurs within S2 stream sur- 
face, three-dimensional flow field must be abominable 
and design progress should be restart. Pressure loss dis- 
tribution could present an objective orientation whether 
the duct wall’s curvature design is acceptable or not. 


Total pressure coefficient Cpt is defined here as: 


In which Pt is local total pressure, Pty the mass aver- 
aged inlet total pressure of ITD respectively, Cpt 
represents total pressure loss. Larger value of Cpt denotes 
that higher total pressure loss will be occurred. 

Figure 8 and Figure 9 are Cpt and Ma contours of dif- 
ferent models, only A8 and A40 are presented. The 
dashed white lines imply the stream-wise planes, S3 sur- 
face, C1 to C5 correspondingly. C1 to C5 planes are ap- 
proximately perpendicular to streamlines along ITD, 
which is intended for analyzing the detailed flow field 
correspondingly. 

Within ITD, ultra high total pressure loss spots mainly 
concentrate on the hub wall while the RA of ITD de- 
creases (red dashed line ellipse areas in Figure 8). This 
means a thicker hub boundary layer on inner wall of 
smaller radial offset ITD (Noting that no separation hap- 
pened within all the models except A45). The inner wall 
curvature effects the low velocity region, which moves 
gradually towards the first bend of ITD in higher ITDs 
(dashed blue ellipses in Figure 9). As increase of ITD’s 
mean radius, the curvature of ITD wall differs. This 
causes variable positions of the peak adverse pressure 
gradient along inner wall. As the result, the main low 
energy fluid congregated near the first bend of ITD. 

The feature of pressure loss near casing is different 
from that near hub. The high total pressure loss area near 
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casing keeps almost the same. It extends from swirl blade 
to LPT nozzle (Figure 8). However, more evident and 
lower velocity area within ITD appears (Figure 9), and 
locates near LPT nozzle for all models. As RA increases, 
the adverse pressure gradient increases along outer wall 
within ITD, and the peak magnitude region keeps nears 
LPT nozzle. 

Boundary layers on duct walls are influenced by RA 
severely according to above discussion. In order to fur- 
ther find out how the shearing stress on walls is affected, 
streamlines and SP coefficient, Cps, are presented simul- 
taneously. 

SP coefficient, Cps, is defined as: 

Cps = (Ps-Pso)/(Pto—Pso) 

where Ps is the local static pressure, Pj) and P.o the 
mass averaged total and static pressure at inlet of ITD 
correspondingly. Cps could distinguish how the SP field 
is changed by various factors while ruling out velocity. 


Fig. 8 Cpt of S2 Stream Surface 


A8 Ma 


Fig.9 Ma of S2 Stream Surface 
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In all the ITDs, LPT nozzle locates before the 2™ bend 
of the ITD, while the 1“ bend starts from right down- 
stream the ITD inlet section (the starting position of the 
figures). From Figure 10(a), it can be easily found that 
completely different limited streamlines appear on the 
hub wall. For A8, smooth surface limited streamlines 
near inlet concentrated near the hub wall, while stream- 
lines before LPT nozzle are nearly perpendicular to the 
axial direction, which means an extremely large circum- 
ferential velocity compared with the axial velocity. Such 
phenomenon indicates that relatively extremely thick & 
low velocity boundary layer concentrated near the ITD’s 
hub wall. 

Variation of the saddle point location, with increasing 
of the RA, near the LPT nozzle’s leading edge, implies 
different inlet angle performance. Comparatively, limited 
streamlines’ skewness near ITD inlet is more overt when 
the RA increases. Span-wise pressure gradient which is 
changed by the 1“ and 2" bend of the ITD induce such 
difference. With larger RA, higher span-wise pressure 


(b) Limited Streamline and Cps on Shroud 


Fig. 10 
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gradient can be found at both the 1“ and 2™ bend within 
the ITD, which causes the low momentum boundary 
layer flow that accumulated near the hub, at the 1“ bend 
position, mixing with the main-stream flow. Such kind of 
flow movement may result in less low momentum flow 
accumulating near the hub wall, and hence reduces the 
streamlines’ skewness. Cps contours reaffirm such varia- 
tion again. 

Similarly, the flow movement close to the outer wall 
differs from each other. The limited streamline on shroud 
of A8 is smoother comparing with its counterpart A40 
case, which is formed by the smaller stream-wise pres- 
sure gradient near the first bend Comparing the limited 
streamlines on hub and shroud, it is influenced more eas- 
ily near inner wall. 

Figure 11 is Cpt of S2 stream surface at 50% span. 
The dissipation of trailing edge shed vortex generated by 
upstream swirl vanes can be observed clearly. However, 
since RA is various in different cases, which means that 
the actual length (not axial length) of the ITD is also 
modified. So even the swirl vane’s design remains un- 
changed in all cases, different circumferential positions 
of upstream trailing edge induced wake can be scruti- 
nized at the inlet of the LPT nozzle. Under its influence, 
incidence angle of LPT nozzle no doubt will change. As 
the result, LPT nozzle’s performance will be definitely 
changed. So it is necessary for designers to optimize the 
blade profile as well as the wall curvature, especially the 
parts near the ITD’s outlet section, to improve down- 
stream turbine’s performance. 

Finally, total pressure losses of all the models (except 
A45) at different planes (Figure 8) are compared by Cpt 
in Figure 12. C1-C2 means pressure loss from Plane C1 
to C2. Apparently, pressure loss is the severest between 
C1 to C2 for all models, near 1“ bend of the ITD. At this 
region, upstream wake is strongly dissipated to increase 
pressure loss. After C2, pressure loss is generated by 
wake dissipation and boundary layer accumulation. 


Fig. 11 Cpt at 50% Passage Height 
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Fig. 12 Cpt of different models 


While it is evident that pressure loss exhibits weak rela- 
tionship relatively between pressure loss and RA. Larger 
RA will not lead severer pressure loss consequently. For 
example, Cpt of A20 is the least than all other models. So 
it seems like that RA is not the only decisive factor for 
pressure loss, while other reason should be considered. In 
this paper, wall curvatures are modified by RA, while 
length and area ratio are constant. As a result, in order to 
decrease loss, it is strongly advised that optimization of 
walls should be taken after RA is decided. 


Conclusions 


In this paper, authors investigate the influence of RA 
on upstream blade and ITD with numerical method. On 
the basis of the same original model coming from a real 
engine, seven ITDs are gotten to be simulated. However, 
separation is found in A45 because of too large RA. So it 
is advised that RA should be smaller than 45 to prevent 
separation. After analysis, several coclusions are drawn. 

Larger RA will lead high loss area (also be the low Ma 
area) near hub to move towards ITD inlet and increase 
the risk of separation within ITD. 

At the inlet section of ITD(outlet of swirl vanes), 
relatively lower Ma area distributes near inner wall 
region and lower yaw angle area seperates near outer 
casing when RA increases to infer a stronger potential 
effect of ITD. So it is issential for designers to modify or 
optimize HPT rotor blade profile while improving ITD. 
More specifically, the outlet angle, as well as stagger 
angle should be larger to keep a relatively constant angle, 
especially near shroud. 

With various RA, the limited streamlines on walls 
present individual characteristics. More specifically, they 
are more twisted near the first bend to indicate the trend 
that the flow field near the first bend is effected most. 
However, the total pressure loss analysis shows that the 
loss is not necessarily larger in ITD with higher outlet 


radius. So an acceptable loss level with large RA is 
possible. Cpt is superiorly high when the ITD is too 
aggressive. 

As a conclusion, the RA’s selection plays a key role in 
making sure a well-designed ITD. Since the flow thread 
through the ITD will decide the downstream flow field’s 
quality and hence the downstream part’s aerodynamic 
characteristic, so the flow movement characteristic 
should be kept in mind when the ITD is studied and 
ready to be used in an engine. An important issue of this 
paper is the total pressure loss could be controlled in 
larger RA ITDs, while optimization is essential. 
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